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NON-LINEAR RESONANCE 
OF FLUIDS 

IN A CRYSTAL GROWTH CAVITY 


INTRODUCTION 


In the microgravity environment, the effect of gravity on fluid motion is much reduced. 
Hence, secondary effects such as vibrations, jitters, surface tension, capillary effects, and electro- 
magnetic forces become the dominant mechanism of fluid convection. Numerous studies have 
been conducted to investigate fluid behavior in microgravity with the ultimate goal of developing 
processes with minimal influence from convection. Industrial applications such as crystal growth 
from solidification of melt and protein growth for pharmatheuticaJ application are just a few 
examples of the vast potential benefit that can be reaped from material processing in space. 

However, a space laboratory is not immune from all undesirable disturbances and it is 
imperative that such disturbances be well understood, quantifiable, and controlled. Non-uniform 
and transient accelerations such as vibrations, g-jitters, and impulsive accelerations exist as a 
result of crew activities, space vehicle maneuvering, and the operations of on-board equipment. 
Measurements conducted on-board a U. S. Spacelab showed the existence of vibrations in the 
frequency range of 1 to 100 Hz. 1> with a dominant mode of 17 Hz and harmonics of 54 Hz. The 
observed vibration is not limited to any coordinate plane but exists in all directions. Similar 
situation exists on-board the Russian MIR space station. Due to the large structure of its design, 
the future International Space Station will have its own characteristic vibration spectrum. 

It is well known that vibration can exert substantial influence on heat and mass transfer 
processes, thus hindering any attempts to achieve a diffusion-limited process. Experiments on 
vibration convection for a liquid-filled enclosure under one-g environment [2,3] showed the 
existence of different flow regimes as vibration frequency and intensity changes. Results showed 
the existence of a resonant frequency, near which the enhancement is the strongest, and the 
existence of a high frequency asymptote. Numerical simulations of vibration convection have been 
conducted by Yurkov^, Fu and Shieh^, and by Wang^ These analyses considered a two- 
dimensional air-filled cell under weightlessness condition and showed results similar to those of 
the experiments. 

It is not yet known whether resonance convection can be triggered by g-jitter alone or 
whether it requires the interaction of g-jitter with other convective forces in low gravity. An order 
of magnitude analysis, however, can be used to show the dependence of the resonance frequency 
on the fluid Prandtl number. Even though the onset of resonance convection may depend on other 
factors, results indicates that fluids with low Prandtl numbers are more susceptible to resonance 
than those with high Prandtl numbers. The current study is aimed at gaining additional insights to 
this problem using germanium as working fluid. Germanium was chosen for this analysis because 
of its common usage in solidification process and its relatively low Prandtl number (Pr = 0.02). 
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analytical MODFT 


The ^viCT-Slokes equations for an incompressible fluid with Boussinesq approximation 
we used for tins analysis. These equations can be non-dimensionalized using fcZSc 
velocity based on thermal diffusion, a, and the dimension of the cell, L. For a two-dtaVnfiond 

thefin" 1 Smuso '? aI v,bration is “‘“"S the same direction as the residual acceleration, 
the following equations cam be used: acceleration. 


o 

II 

(1) 

cP 

* + ^ +F ^F = -ir prvc/ 

(2) 

^- + U— + V— = -^- + Pr v 2 F 
& ax dY cY 

+Pr[^?<ar + G(o sin (zeo)]0 

(3) 

St ax ar ax 2 ar 1 

(4) 


Here U and V are the non-dimensional velocity components, Pr the Prandtl Number Ra the 
^yleigh number, G the Grashof Number, and © the non-dimensional circular frequency Notice 
at the forcing term m equation (3) is proportional to the product of Prandtl number the 
emperature differential, and the buoyancy effect of acceleration. The buoyancy term consists in 
general, of two components, one due to the residual acceleration and one due to the oscillation 
o°.dlt J he C ° nStant accelerat,on IS expressed in terms of the Rayleigh number, and the applied 

of the product of the the ^ ^ - 


NUMERICAL INVESTIGATION 


e- •. ?J uations O-f) Were solved numerical| y us “>g the NASA developed two-dimensional 

which 1 enCe Nav l er i' oi ; es solver (FDNS2D) [Ref 7], The working fluid is germanium 
whtch ,s a representauve flutd used in solidification from melts. The properties of germamum 
along with that of water and air are given in Table 1. germanium 


Property 


Air Germanium 


Kinematic Viscosity (m2/s) 
Density (kg/m2) 

Thermal Diffusivity (m2/s) 
Expansion Coefficient (1/K) 
Prandtl number 


9.8260E-07 

9.9740E+02 

1.4490E-07 

2.7720E-04 

7.00 


1.6000E-05 

1.1650E+00 

2.2860E-05 

3.3000E-03 

0.70 


3.3160E-07 

6.1000E+03 

1.5120E-05 

1.2600E-04 

0.02 


Table 1. Summary of Fluid Properties 
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The calculation domain of the cavity is divided into a finite difference mesh of 3 1 by 3 1 . 
The fluid was initially assumed in a quasi-static condition with a constant temperature difference 
of 100k maintained across two opposing surfaces. At time zero, a continuous vibration 


acceleration of constant amplitude is applied in the direction normal to the temperature gradient. 
Transient solutions were obtained for different vibiation frequencies using a time marching 



L 


W©4 




M i flixrh C. 


; : ; • - • • 
iii : ; : ' 

■ < i « 1 



‘SiJe. 

i : ; : ; 

; ; ! : 1 ! 

1 J t 1 1 1 * 

1 ! , 1 1 1 

T c 

’ 1 i 1 ' 

l , 1 , *• 1 

. ; . j 1 « | 

— 1 > 





Med ri«*. 


Figure 1 . Geometry of Cavity 


Figure 2. Average Nu Plot for Germanium 


RESULTS & DISCUSSIONS 


Numerical results for germanium (Pr = 0.02) in a cavity of 1 cm squared in a zero gravity 
environment will be reported. A Grashof number of 10 + ^ was used in the analysis. Figure 2 shows 
the existence of several different flow regimes as the vibration frequency increases. Expressed in 
terms of the averaged Nusselt number, the result indicated an increase of heat transfer from the 
hot wall to the cold wall at low vibration frequencies. This heat transfer rate peaks out near a 
resonance frequency, then went through a transition regime of instability. Finally it settles to an 
asymptote at very high frequencies. This trend is very similar to that of air, shown in Figure 3 
from a previous study [Ref. 6], However, the transition between the flow regimes occurs at much 
lower frequencies, and the enhancement of heat transfer rate is less prominent for germanium. 

The resonance regime for germanium is found between the non-dimensional frequencies of 3 and 
10, which translates to a frequency of 0.247 Hz for a cell of 1 cm. This frequency is much smaller than that 
of air which exists in the range of 10 to 30 Hz. This difference is due to the difference in fluid 
properties, especially due to the fluid viscosity and fluid Prandtl number. A plot of the steady state 
temperature at a point near the top of the mid-plane is shown in Figure 4. It shows that at lower 
vibration frequencies the temperature variation is high, and its mean value is close to the medium 
value between the two side walls. As the frequency increases, the temperature variation becomes 
smaller and the mean value moves toward a lower value. A review of temperature data at several 
1 ( 


to the direction of vibration, with the top boundary always at higher temperature than the lower 
boundary. This has been attributed to the persistence of flow as a result of the initial impulse. 

Time variation of the averaged Nusselt number and temperature at three points along the 
mid-plane for a few selected frequencies are given in Figures 5 and 6, respectively. The three 
frequencies selected co = 3, 10, and 100 represent regimes of transition, resonance, and high 
frequency vibration convection. Notice the existence of the double peaks in Nusselt number 
during each cycle of vibration for co = 3 in Figure 5. This is due to the complete reversal of the 
temperature about the mid-plane of the cavity during an oscillating cycle of vibration, as shown 
clearly in the top portion of Figure 6. 
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Figure 3 . Average Nu Plot for Air 
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Figure 5. Variation of Nusselt Numbers for 
co = 3 (top), 10 (middle), and 100 (bottom) 



Figure 4. Temperature Plot for Mid-Plane 




Figure 6. Variation of Temperature for 
co = 3 (top), 10 (middle), and 100 (bottom) 


LIV-4 






CONCLUSIONS & RECOMMENDATIONS 


Results of a numerical simulation of vibration convection of a crystal growth cavity in zero 
gravity showed the existence of different flow regimes at different vibration frequencies. Results 
showed that the division of the flow regimes and the peak heat transfer rate dependent on the 
Prandtl number of the working fluid. A complete reversal of temperature profiles about the cavity 
mid-plane at low vibration frequencies leads to heat transfer rate which oscillates at twice the 
vibration frequency. At higher vibration frequencies, the fluid near the top portion of the cavity is 
consistently hotter than the fluid near the bottom, resulting in an internal heat flow in the direction 
normal to the applied temperature gradient. At very high frequencies, the flow reaches an 
asymptote, and the heat transfer rate becomes independent of the vibration frequency. 

Additional analyses are needed to investigate the effects of cavity geometry, e.g. 
cylindrical shape, or one with different length to height ratios. Effects due to non-zero gravity also 
need to be considered. Experimental programs aimed at verifying the numerical results are under 
consideration. A flight experiment will eventually be needed to verify the results of vibration 
effects on fluid convection in micro-gravity conditions. 
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